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ABSTRACT 



A digital data acquisition system was implemented for the study and collection 
of scattered waveforms from various simple objects. Reflected waveforms were 
extracted from interference with other signals by a simple subtraction technique. 
Backscattered and bistatic reflections were observed and recorded. Various effects 
on scattered waveforms, including target composition, pulse length, and frequency, 
were experimentally observed. Bistatic scattering results from a stryofoam sphere 
at a ka value of 10 were in agreement with theoretical predictions for a rigid sphere. 
Bistatically measured data at 70 kHz from a styrofoam cylinder of length 20 cm and 
diameter 2.5 cm followed a trend for bistatic angles up to 150 degrees but became 
amorphous as the bistatic angle approached 180 degrees. 
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I. TNTRODIJCTTON 



Any object with a characteristic impedance different from the surrounding 
fluid will affect an acoustic wave impinging on it. The wave will be transmitted 
through the object and scattered from it. Flexural waves, compressional waves, and 
creeping waves are examples of some of the possible scattering and reflection 
mechanisms [Ref. l:p. 301]. The "characteristic impedance" equals the density of a 
material multiplied by the speed of sound in the material. It is one measure used to 
calculate reflection factors for pressure waves. Other properties that affect the 
scattering properties of an object include frequency, shape of the object, and pulse 
length. 

Ever since the development of the submarine as an important force in naval 
warfare, much effort has been directed to improving their detection. The same 
holds true for the submarine detecting of surface ships or other submarines. The 
knowledge required for detection include knowing the source level of the target, 
propagation loss, noise level , and lastly target strength, the subject of this paper. 
These parameters are used in the well-known sonar equations [Ref. 2;p. 411]. 
Target strength, TS, as used in the sonar equations is defined as: 

TS=101ogil (1) 

1q 

where If is the reflected intensity one meter from the acoustic center of the target 
and lo is the incident intensity at the target without the presence of the target. This 
equation is general, giving no stipulation as to the direction of the reflected 
intensity. Typically, signals in active sonar are sent and recevied by the same 
platform, thus the most common geometry and , therefore, the most studied, is the 
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backscattered monostatic case where the reflected intensity is measured in the 
direction back toward the source. Measurement along any other direction 
constitutes a bistatic geometry. Relatively little effort has been given to the study of 
bistatic target strength until recently when the bistatic geometry was recognized as a 
potential method for improving underwater detection of submerged targets. 
Figures 1 and 2 illustrate the monostatic and bistatic geometries. 



Source 

and 

Receiver 




Target 



Figure 1. Backscattered Monostatic Geometry 
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Receiver 



Figure 2. Bistatic Geometry 
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Some work has been done in the radar community on estimating the bistatic 
radar cross-section of objects. Crispin and Siegel [Ref. 3:p. 156] present a theorem 
for the approximation of radar bistatic cross-section: 

In the limit of vanishing wavelength, the bistatic cross-section for 
transmitter direction tc and receiver n is equ^ to the monostatic cross-section 
for the transmitter-receiver direction k -i- n with I: n for bodies that are 

sufficiently smooth. 



(Figure 3 illustrates the geometry as stated in the theorem.) 




Figure 3. Radar Bistatic Cross-Section Theorem Geometry 

Their bistatic theorem has its limitations wHen applied to complex shapes; The 

bistatic cross-section can be greater or less than the monostatic case due to 

/■ 

structurally-related interferences. Another limitation is the size of the target with 
respect to the incident wavelength. Spheres were analyzed using various ka values 
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where where a is the radius of the sphere and k is the wave number. It was found 
that the theorem held more closely for higher ka values, but in all cases the theorem 
sharply deviated from experiment as the bistatic angle approached 180 degrees. As 
an illustrative example, they presented data on a jet aircraft and a ballistic-missile 
shape [Ref. 3:p. 174-177]. Data were in fairly good agreement with their theorem 
as long as the bistatic angle did not get too large. 

In transition from the radar community to underwater acoustics, Urick [Ref. 
l:p. 302] implied that the bistatic radar cross-section theorem might apply in a 
underwater condition, in which case the bistatic target strength was equal to the 
monostatic target strength taken at the bisector of the bistatic angle. Urick 
emphasized there was a lack of data to prove or disprove such a theory. It would be 
very convenient if such a theorem existed since there is plentiful target strength data 
for monostatic geometries. Whether or not the high frequency limitations of such a 
theorem would preclude its tactical use is another question. 

The interest in bistatic target strength has been generated by the possibility of 
increasing detection of submerged and surface units using this technique. For 
example, one can easily visualize having a surface ship using its active sonar while a 
submerged submarine, using passive sonar, listens for reflections off enemy forces. 

There are several advantages and disadvantages to bistatic employment. Several 
ideas now discussed were taken from Harvey [Ref. 4] who gave a thorough 
discussion on using bistatic employment as a tactical tool in naval warfare. The 
major advantages include the ability to maintain a track of the enemy force using 
passive sonar and the increased amount of data available to detect and classify 
targets. All units besides the active unit can be passive, thereby hampering own- 
ship detection. Detecting bistatic reflections, as well as backscattered reflections. 
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from a target increases the data available for detection and classification of the 
target. The greatest difficulty using a bistatic geometry is the high positional 
accuracy required between the pinger and the listener. This requires prior planning 
for use with surface ships and submarines since they do not normally have 
continuous open communication with each other. Another difficulty is the 
distinguishing of backscattered reflections from bistatic reflections. It may turn out 
that the advantages do not outweigh the disadvantages and the whole idea will prove 
infeasible. 
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IL TARGET STRENGTH 



A. INTRODUCTION 

This chapter presents some definitions and background knowledge essential in 
preparation for experimentation. Some concepts are rather simple but are included 
for completeness. 

B. DEFINITION 

The definition of target strength has already been given, equation (1). 
However, there are several aspects of the definition that warrant clarification. 
First, all measurements are taken in the far field of the target. The far field is 
where the pressure reflected from an object exhibits a 1/r relationship [Ref. 2:p. 
178]. For example, pressure measurements taken one meter from the hull of a large 
ship would certainly be in the near field, so intensity is measured in the far field at 
ranges of hundreds of meters and then extrapolated back to one meter for use in the 
target strength calculations. 

The next question concerns the defintion of the acoustic center. Very rarely are 
we able to know the location of a target’s acoustic center. The acoustic center is the 
point from which sound appears to radiate. It may be anywhere around the target, 
and may lie outside the shape itself. It may also change with the aspect of the target 
relative to the receiver. Exact location of the target’s acoustic center is ignored in 
field measurements, and the geometrical center is used. When measuring pressure 
at hundreds of meters, a few meters error in determination accounted for by the 
acoustic center will not produce a large discrepancy. 
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C. SPHERICAL WAVES 



For a distant source, to a good approximation, the wave will appear locally as a 
plane wave as it encounters the target. Consider an omnidirectional point source: 
Spherically diverging waves are emitted and the pressure amplitude and phase is 
dependent only on the distance r from the acoustic center of the source. An 
outwardly travelling spherical wave can be represented by 




j (wt - kr) 



( 2 ) 



A = pressure amplitude at one meter 
r = distance from point source 
w = angular frequency of signal 
k = wave propagation number 

When the radius of curvature of the wave front is large compared to the biggest 
dimension of the target, the spherical wave can be treated as a plane wave. 

Plane waves are assurred when the dimensions of a target are small compared to 
the source-target distance. One of the goals of small-scale experiments is the 
application of the results to oceanic situations, therefore there is a concern about the 
approximation of plane waves ensonifying a target in a small laboratory tank. 

The geometrical constraints of the water tank limits target-source-receiver 
distances. With small distances between source and receiver and big targets, a plane 
wave may not be a good approximation. Another factor that influences the plane 
wave approximation is the directionality of the source. It is difficult to achieve 
omnidirectional behavior from a finite source and it is impossible if the frequency 
is high. Beam patterns are formed with energy emitted in preferred directions. If 
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the target subtends an angle greater than the beamwidth, the pressure amplitude will 
not be the same for all points on the target. 

Plane waves also have properties that make target strength calculations simple 
to handle. The intensity can be equated to the root-mean-square pressure amplitude 
squared divided by the characteristic impedance [Ref. 2:p. 111]. Since hydrophones 
measure pressure via a sensitivity factor, the output voltage they produce can be 
directly related to intensity. 

D. WAVEFORM EFFECTS 

There are several factors to consider in the measurement and collection of 
reflected waveforms from objects. Some are listed below and each will cause the 
reflected waveform to differ from the incident waveform. 

- target composition 

- target shape 

- pulse length 

- resonance effects 

For the reflected waveform to be identical of the incident waveform, a planar 
rigid boundary must be encountered. For more complicated targets, some of the 
incident energy will be transmitted into the body and eventually reradiated and 
different areas of the target may reflect in such a way to cause mutual interference. 
Many objects have resonances, at which the target will store some of the incident 
energy and reradiate it at a later time causing interference effects. Theoretical and 
experimental results show that for short pulse lengths, reflections are nearly 
independent of target composition and are almost identical to the incident pulse 
[Ref. 5,6,7] while for longer pulses the reflected waveform changes drastically 
depending on the composition and shape of the target. 
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E. EXPERIMENTAL MEASUREMENT 

While conducting the literature search, it was interesting to note how the 
pressure amplitude of the reflected wave was measured. As stated before, 
hydrophones give a voltage that is representative of pressure and this can be equated 
to intensity. The problem lies in the fact that the reflected waveform, in most cases, 
is not identical to that of the incident signal. The experimenter has to specify what is 
measured and what kind of pulse is being used. Response can be catagorized as 
transient or steady-state. Steady-state implies that the whole target has been 
ensonified by the incident pulse allowing for all portions of the target to contribute 
to the reflection and for any resonances to become fully excited. In this case, the 
return will be a combination of all possible reflection paths and modes. A transient 
response simply is not in the steady-state. To give an idea of measurement used, the 
following are quoted from two published articles: 

...five cycles of the sine wave located at a position in the pulse judged to 
be a true representation of the steady state. The average of the ten peak values 
so obtained gave the peak pressure amplitude [Ref. 7] 

After amplification, the amplitude of the back- scattered signal is 
measured with a peak detector [Ref. 8]. 

The main point is that for reflected waveforms, not similiar in shape to the 
incident waveform, some subjective determination is required. Estimation of the 
steady-state must be made or some definite parameter has to be selected such as peak 
amplitude or total energy. In any case, an unique way to convert the amplitude 
measurements to an intensity level is not achievable. Add to this the fact that the 
reflected waveform will be different for different target orientations and 
frequencies, and the problem of translating measurements into an intensity level to 
calculate a target strength becomes comphcated. 
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There presently exists a Target Strength Standardization Program [Ref. 9] 
which was assigned the task of finding some standard measure of target strength 
independent of signal type. The general consensus from the program was to use the 
"Average Integrated Target Strength" which is the average time-integrated energy. 
Furthermore, an Echo Metrics Subcommittee [Ref. 9] was formed to find signal- 
processor dependent metrics for the measurement of target strength. Theoretically, 
the incident energy encountered by the object would be reflected or absorbed and 
eventually reradiated at some later time providing no energy was stored or lost in 
the object. Therefore, it would seem that the average of the integrated energy 
would be an independent metric for all types of signals providing the integration 
time is extended to include all the return signal. Some standardization of target 
strength is being investigated so that data from independent experiments can be 
compared and integrated. 
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III. DATA ACQUISITION SYSTEM 



A. BACKGROUND AND CONCEPT 

The experiments were conducted in a water tank located in Room 019 in 
Spanagel Hall at the Naval Postgraduate School. The tank was approximately 1 .5 m 
wide, 8 m long, and 2 m deep and lined with anechoic rubber material. For the 
targets used, wall reflections were of the same magnitude as the target reflections. 
Usually, it was easier to use pulsed signals and time discriminate for the target 
reflection. However, the small dimensions of the tank did not provide the 
environment to do this effectively for all geometries. To overcome this difficulty, a 
subtraction technique was used employing a digital storage device. The basic 
concept was to subtract the signal received with the target immersed in the tank 
from the signal received with the target absent. The procedure was very simple and 
allowed for any target-source-receiver geometry. 

B. DEVELOPMENT 

A considerable amount of time went into the selection of the digital sampling 
and storage device, the source, and the receiver. The first question that had to be 
addressed was the range of frequencies to be used. The frequency range of 100-200 
kHz was initially selected because it gave wavelength-to-target dimension ratios that 
match oceanic scenarios. The critical specifications for the digital storage device 
are the sampling rate, dynamic range and jitter. Dynamic range was not a problem 
since most digital oscilloscopes have the same voltage range as a typical non-digital 
oscilloscope. The sampling rate must be fast enough to provide a good 
representation of the signal, and assuming a sinusoidal signal about ten samples per 
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cycle would give a good trace. For example, for a 100 kHz signal the sampling rate 
would have to be 1 MHz. Another major concern was the amount of jitter in the 
sampling device. Jitter describes the fluctuation in the triggering times relative to 
the sampled signal. For a system using a subtraction process, it is perhaps the most 
critical specification because there has to exist a common time reference between 
collected waveforms to get accurate results. 

The selection of the digital oscilloscope was based on the availability of 
equipment. Two digital oscilloscopes were evaluated. The Nicolet 3091 was 
analyzed first. Its fastest sampling rate was 1 MHz and had a storage capability of 
4000 points. The Nicolet 3091 had a RS-232 serial interface and was easy to use, 
however it was not acceptable due to excessive jitter at the frequencies of interest. 
The jitter was easily seen by applying a pulsed signal and observing the traces on the 
scope. The digital oscilloscope finally used was the Tektronix 2430 (Tek 2430). It 
is a very capable oscilloscope and has many nice features making it suitable for data 
acquisition. The most important feature that it had was an averaging mode 
selectable from 2 to 256 waveforms. This helped eliminate random noise and 
minimized any fluctuations the electronic equipment possessed. The amount of 
jitter was dependent of the time scale selected. The fastest sampling rate of the Tek 
2430 was 100 MHz which could be used to sample a 10 MHz signal using ten points 
per cycle. The jitter specification given in the operating manual, for the time- 
delayed mode, was .08 x (time scale). The jitter was observed and was significant 
for the time scale desired, however the saving factor was the time averaging mode 
which ultimately averaged out the jitter and provided a sufficiently constant time 
reference. It was very easy to observe this on the oscilloscope because of a 
triggering marker on the signal. Jitter becomes less of a problem for any digital 
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oscilloscope when the number of samples per cycle of signal increases. However, 
the problem then arises from the limited amount of storage in the oscilloscope. A 
smaller number of cycles can be stored when using a faster sampling rate. The Tex 
2430 has a storage capabiltiy of 1024 points, much less that the Nicolet 3091; 
however it also has a time-delay feature which allowed for selective data storage 
(The Nicolet 3091 did not have any time-delay feature). Other features of the Tek 
2430 include cursor control for voltage and time measurements, simple output to a 
Hewlett Packard Think-Jet printer or pen plotter, and a GPIB interface. It is a very 
versatile piece of equipment and was finally selected for use in the data acquisition 
system. 

The selection of the source and receiving hydrophone was also limited to in- 
house equipment. For the frequencies desired, the USRD type F41 transducer was 
the appropriate choice with a nominal frequency range of 15-150 kHz [Ref. 10:p. 5]. 
It is a piston type piezeoelectric transducer with a face diameter of 10 cm. The 
choice of the receiver was a Celesco LC-10 hydrophone. It has a wide nominal 
frequency band, is small and light making it easy to move around, and has good 
sensitivity. 

C. SYSTEM SET-UP 

1. Electronic Equipment 

The F41 source had to be driven near 200 V rms for target reflections to be 
seen. An Instruments, Inc. Model L2 power amplifier was used to drive the source 
and its selection was again based on what was available in the laboratory. It allowed 
for impedance matching and has a voltage-divided 10 mVA^ port to monitor output 
voltage. The output of the receiver was passed through an Ithaco Model 1201 Low 
Noise preamplifier which has bandwidth setting capability. This was very 
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beneficial because there was a lot of 60 Hz noise which had to be eliminated by 
bandwidth filtering. The signal out of the preamplifier was sent to a Kikusui 
COS5060A oscilloscope and the Tektronix 2430 digital oscilloscope. Having a non- 
digital oscilloscope made observing and measurement of the waveforms easier to 
perform on the Tek 2430. A Hewlett Packard 3314A function generator was used 
to provide the source signal. It also provided the triggering signal which was sent to 
the external channels of the oscilloscopes. In the pulse mode, the HP 3314A sends 
out a 3 V transistor-transistor logic signal via its trigger port. The pulse repetition 
rate of 10 ms was sufficient to allow reflections to diminish before the next pulse. A 
Hewlett Packard 9000 Series 300 desk top computer was used as the data collection 
and processing unit. It was connected to the Tek 2430 via the HPIB bus. The 
transfer of waveforms stored on the Tek 2430 was controlled by software 
programs. Data was ASCII encoded from the Tek 2430, transferred, decoded into 
actual voltages, and then stored onto disks in binary data form. A schematic 
diagram of the electronic system is shown in Figure 4. 

2. Physical Set-up 

To measure all possible bistatic angles, the scattering pattern would have to 
obtained in three dimensions. This was not done in this work. Bistatic data were 
limited to the horizontal plane determined by the source, receiver, and target, with 
the axis of the target in the same plane. 

The source, the F41 transducer, was mounted on a 3/4-inch aluminum rod 
attached to a rectangular support plate. The rod was fixed to the plate by means of a 
friction fitting secured by a hose clamp. Target, source, and receiver were 
submerged one meter below the surface of the water and the rods which supported 
them were marked for this depth. The source with its support rested on another 
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Figure 4. Data Acquisition System 




support built from aluminum L-beams that spanned the width of the tank. Lead 
bricks were used to hold the source once it was set in place and a level was used to 
ensure the source face was vertical. The receiver was mounted onto a sliding track 
specially built for the tanks. Adjusting screws were added for ease of positioning 
the receiver vertically. The initial intent was to design a rotating boom to hold the 
receiver and support the target from its center. A device was constructed, but too 
much receiver movement occurred between the loading and unloading of targets. 
Though the boom was never used, the device was still used for the suspension of 
targets and had a mechanical angle indicator on its top for angle measurement. 
Spherical targets were suspended by a thin cable and the cylindrical target by a 3/8- 
inch rod. 

The direction of the source was set by aligning the source, receiver, and 
target along a straight line with the receiver between the source and target. Source, 
receiver, and target depth were one meter. Vertical alignment was ensured by 
keeping the source and receiver vertical by use of a level. The frequency was set to 
200 kHz to establish a narrow major beam and the source was rotated until 
maximum signal was seen by the receiver. This ensured the target was being totally 
ensonified by the major lobe of the source. 

Once the source direction was set, the receiver was moved to the desired 
position. The angle between the source, target, and receiver was determined by 
setting the receiver’s track perpendicular to the source-target axis and then using 
trigonometry relationships. After the desired geometry was set, care was taken to 
ensure the target, source, and receiver were stationary. It was critical that all 
elements be stationary when waveforms were collected. 
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3. Data Collection 



The method of data collection and storage followed the sequence: 

- set up desired source-receiver-target geometry 

- warmup all electronic equipment >1 hr 

- ensure connections as illustrated in Figure 5 

- load ”Tex" onto computer 

- set up the Textronix 2430 oscilloscope 

- B trigger mode to "runs after" 

- A trigger 

- mode; normal 

- coupling: AC 

- source: external 1 

- slope: positive 

- level: 800 mV 

- trigger position: 1/8 

- data encoding: ASCII 

- GPIB address: 4 

- vertical mode: channel 2 

- select time delay to obtain desired time window 

- set acquisition mode to 256 averaging 

- adjust channel 2 voltage scale/Ithaco gain appropriately 

for the signal 

- for each new signal, press the acquire button 

- wait sufficient time for averaging to take place 

- press the save button 

- run the "Tex" program 

The "Tex" program controlled the transfer and storage of data from the 
digital oscilloscope. Inputs from the user were file names and file descriptions. 
The transfer process, including writing onto the disk, took approximately 10 s per 
waveform. A listing of the program is in Appendix A. 
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The computer program "Plot" was used to evaluate the stored data. It 
took two files and subtracted them, and then plotted the results on the CRT 
terminal. Hard copy output could be obtained by porting to a pen plotter or Think- 
Jet printer. Input from the user required file names and where printout was to 
occur. The plotting routine was written so that the plot’s parameters were easily 
changeable. A listing is enclosed in Appendix B. 

The method of data collection was a slow and somewhat tedious process 
since all movement of target and receiver ware manually controlled. Data 
readings were also recorded manually. 

4. Validation 

The question of whether the data acquisition system could correctly 
extract a target reflection immersed within other signals had to be answered. The 
very simple procedure of subtracting background noise works only if there is a 
common time reference for "target in" and "target out" waveforms. For ease of 
discussion, "background noise" is used for any signal, random or coherent, that is 
not the desired signal and corresponds to the "target out" condition. 

The procedure to validate the system began with obtaining a geometry 
which allowed a direct pulse and reflected pulse to be completely separated by 
time discrimation with no interfering wall reflections. Figure 5 is such an example. 
The larger pulse was the direct path pulse and the second pulse is the pulse 
reflected from a styrofoam sphere. The frequency was 66.5 kHz and the pulse 
length was 20 cycles. The reflected pulse is shown enlarged in Figure 6. The 
peak-to-peak voltage of the reflected pulse was measured to be 68 ± 2 mV using 
the cursor controls of the Tek 2430. 
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The pulse length was then increased to 120 cycles so the direct pulse would 
be interfering with the reflected pulse. No other changes were made. The 
interference of the two signals is shown in Figure 7 and was stored in file A. Next, 
the sphere was removed and the background noise sampled and recorded in file B. 
The background noise in this case was mostly the direct pulse. The last step was to 
subtract file B from file A to get the reflected signal from the sphere. The result is 
shown in figure 8. The waveform begins at the same time in comparison with 
Figure 6 and the only difference is that the pulse length was longer as expected and 
the curve slightly more choppy in form. Measured voltage was 73 + 4 mV. The 
larger uncertainty arose due to the combination of small jitter and amplitude 
fluctuations of the HP 3314A. The uncertainty was estimated beforehand by taking 
identical signals directly from the function generator to the Tek 2430, recording 
them in separate files, and subtracting them and observing the results. It was found 
that peak-to-peak voltage uncertainty was approximately four digitizing levels 
regardless of the voltage scale used (a digitizing level equals the voltage per division 
divided by 50 for the Tek 2430). Given that this error would always be present, 
there was good agreement between Figure 6 and Figure 8, and the method of 
obtaining data was considered validated. 

As with any signal processing system, a high signal-to-noise ratio is 
desired. For any digital subtraction routine, the desired signal has to be large 
enough relative to the noise so it will not disappear in the digitizing level of the 
sampling device. For example, if noise plus the signal were 10 volts using 512 
digitizing levels or bins, each bin represents 19.5 mV. If the signal were small, i.e., 
1 mV, the signal level would be much smaller that the digitizing level and would be 
lost in the subtraction routine. 
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Figure 5. Direct and Reflected Pulse 
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Figure 6. Reflected Pulse 
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Figure 7. Interference Between Direct and Reflected Pulse 
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Figure 8. Reflected Pulse, N=120 
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Keeping the source and receiver positions fixed relative to each other was 
critical. Small movement between the two greatly affected the phase of the received 
signal and had a large effect on the results. This effect was verified in preliminary 
work. The initial intent was to place the receiver and support the target on the same 
platform. However, the removal of the target for the background noise 
measurement caused flexing in the platform and produced source-to-receiver 
movement giving incorrect results in the subtraction routine. Figure 9 illustrates 
the effect of source-receiver movement on "target-in" and "target-out" waveforms. 
The solid trace represents the target-in waveform and the dashed line represented 
the target-out waveform. The reflection began at 35 usee. The dashed and dotted 
line is the the result from the subtraction of the two waveforms and it was clear that 
in this case, it did not give the correct results. The region between 0 and 35 usee 
should have cancelled to zero. The non-zero result indicates a phase change caused 
by relative source-receiver movement. The phase change subsequently affected the 
amplitude of the reflected signal from 35 to 100 usee. It was easy to see that any 
phase change brought about by small source-receiver movement would cause 
erroneous results. To eliminate this problem, a separate platform was subsequently 
built and used for the receiver. 

The output of the power amplifier was also checked in preliminary work 
for fluctuations. An arbitrary source-receiver geometry was selected without any 
target in the tank. Several waveforms of the direct pulse were stored at various 
times after initial start up. Amplitude fluctuations existed shorty after start-up. 
Ideally, successive waveforms subtracted from each other would cancel out to give a 
zero response. If amplitude fluctuations existed a false signal would occur and this 
in fact did occur after five minutes of warm-up as Figure 10 illustrates. Another 
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figure 9. Relative Source-Receiver Movement Effect 



test taken at 15 minutes after stan up produced better but still unsatisfactory results. 
After one hour of warm-up the amplitude fluctuations subsided and a zero respose 
was obtained. 




Figure 10. Waveform Amplitude Fluctuation Effect After 10 Minutes 

of Warm-up 
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IV. FAR FIELD DETERMINATION 



The measurement of the far field of several objects was conducted. An equation 
that gives the approximate beginning of the far field is given by [Ref. 2:p. 188]: 







( 3 ) 



L = the longest dimension of object 
X = wavelength 

It is always best to determine the far field by measuring pressure amplitude as a 
function of distance from the object. The equation for a diverging spherical wave, 
as described by equation 2, can be manipulated by taking the log of both sides to 
transform it into a linear equation with a slope of minus one. The far field can be 
determined experimentally when the observed slope of the curve starts to deviate 
from minus one. 

The far field of the F41 was determined to begin at approximately 100 cm and 
70 cm for 176 kHz and 70 kHz, respectively. The far field for the sound reflected 
from the styrofoam cylinder positioned broadside to the source was determined to 
begin at 40 cm for 70 kHz. The three-inch diameter styrofoam sphere’s far field 
began at less than 20 cm at 70 kHz. Figures 1 1 through 14 are graphs of the data 
used to determine their far fields. Lines of slope equal to minus one are drawn on 
the graphs. 
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Figure 11. Receiver Voltage vs Distance, F41, 176 kHz 
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Figure 12. Receiver Voltage vs Distance, F41, 70 kllz 
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Figure 13. Receiver VoUage vs Distance, Cylinder, 70 kHz 
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Figure 14. Receiver Voltage vs Distance, 3-Inch Dianieter Sphere, 70 kllz 



V. SPHERES 



A. INTRODUCTION 

To achieve an understanding of the complexities involved in analyzing 
scattering effects of various objects, it was instructive to begin with a very simple 
and symmetrical object, a sphere. Several solid metal spheres were used. Spheres 
were aluminum and brass, but, the exact composition was not known or investigated 
since they were to serve only as preliminary tools in studying the general properties 
of waveform scattering. 

Since metal spheres are not perfectly rigid, they are capable of absorbing 
energy and subsequently reradiating it. Also, they possessed resonances. Many 
experiments in the past 30 years have been conducted with spheres and long 
cylinders, and the Resonance Scattering Theorem (RST) has evolved from this work 
[Ref. 5,6,7,8,11,12,13,14]. This theorem predicts scattering patterns which entail 
both the rigid-like properties and the resonance properties of an object. It was not 
the intent to examine this theorem in detail but only to observe some of its 
predictions experimentally. 

B. MONOSTATIC FORM FUNCTION 

The goal of replicating some of the data in the literature was sought to establish 
confidence in the experimentation process. Neubauer, Vogt, and Dragonette [Ref. 
7] conducted an analysis on elastic spheres using steady-state signals and observed 
the backscattered signal form function from an aluminum sphere as a function of ka. 
The reflection form function is defined as: 
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Ps = scattered pressure amplitude measured at r 
po = incident pressure amplitude at sphere center 
a = sphere radius 

Their experimental data were plotted and compared against theoretical results 
and good agreement occurred. 

Part of Neubauer's experiment was repeated. A four-inch diameter aluminum 
sphere was ensonified by 50 cycle pulses from 100 to 130 kHz, and the backscattered 
pressure amplitude was measured. Time discrimination of the reflection was used. 
The amplitude of the reflected pulse was taken at what was judged to be the steady- 
state part of the signal. The receiver was 1 .60 m from the center of the sphere, 
placing the receiver well within the far field of the sphere. 

The wavelength used in the calcuation of ka was dependent on the speed of 
sound in the water, which in turn was affected by temperature, pressure, and 
salinity. For the fresh water tanks, salinity was ignored. The speed of sound was 
based on an empirical formula that applied to distilled water [Ref. l:p. 107]. 
Although the water was normal tap water, the distilled water formula was a close 
enough approximation [Ref. 15]. Water temperature of the tank was a consistent 17 
degrees Celsius throughout the entire experimentation phase. At this temperature 
and normal atmospheric pressure, the speed of sound was calculated to be 1473 m/s. 
Experimentally, two LC-lO's were placed on a track 180 cm apart and the 
successive arrival times of the pulse were measured via the Tek 2430 with the 
cursor controls. The speed of sound thus was measured to be 1480 ± 20 m/s, placing 
1473 m/s well within the experimental uncertainty, and this value was used as the 
speed of sound in water. 
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The shadowing effect of the LC-10 in front of the F41 source was taken into 
account by placing a dummy LC-10 at that position and recording incident pressure 
at the sphere center location. Additionally, since the LC-10 was not perfectly 
omnidirectional, the support rod was marked to maintain a constant receiver aspect 
with the target. Placing the LC-10 in front of the F41 also caused additional 
reflections to occur. The direct pulse was reflected by the LC-10 back to the 
source's face and again reflected there back to the LC-10. More "noise" was present 
when measuring the monostatic case but appropriate location of the receiver 
eliminated any interference from those reflections. 

The data matched Neubauer's results well, and the form function as a function 
of ka is plotted in Figure 15. Data from Neubauer's experiment is also plotted in 
Figure 15. The only difference between results was a small shift in ka of 
approximately 0.2-0. 3, not surprising since the compostion of the sphere was 
unknown. The minima corresponds to resonant frequencies of the sphere, where 
the specularly reflected wave is interfering with the resonantly radiated wave from 
within the body. 

C. WAVEFORM EFFECTS 

The effect of a sphere's resonance was examined in more detail. A four-inch 
diamter aluminum sphere was ensonified by a 30 cycle pulse at, above, and below a 
resonant frequency of 82.7 kHz. At resonance, the ringing of the sphere is clearly 
seen in Figure 16 when the direct pulse ends. Pulse length was approximately 367 
|is long. The sphere's stored energy is simply being reradiated. Figures 17 and 18 
show the waveform response 2 kHz above and below resonance. They are 
decisively different and clearly illustrate the resonance effect. 
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Figure 15. Backscattcred Form Function vs ka, Aluminum Sphere 



The effect of sphere composition on reflected waveforms was observed. A 
four-inch diameter brass sphere was ensonified by the same frequencies as the 
aluminum sphere. Brass has a characteristic impedance approximately 2.5 times 
that of aluminum. Figures 19, 20, and 21 show the reflected waveforms having 
frequency dependence that are markedly different from the aluminum sphere. 

The effect of pulse length was observed using a 4-inch aluminum sphere at and 
off a resonant frequency. Figure 22, at resonance, and Figure 23, off resonance, 
show reflected wavforms using a short pulse of 2 cycles. They appear very similar 
in shape. Figures 24 and 25 show the results for a pulse length of 5 cycles. Again, 
the results of the two are similar. When pulse length reaches 15 cycles, the effects 
of the resonance begins to be clearly observed, as Figure 26 and 27 illustrate. 

For comparison, the same procedure above was performed on a four-inch 
diameter brass sphere. Results are illustrated in Figures 28 through 33. The brass 
sphere did not show any resonance effect as determined in the preceding 
composition experiment for these frequencies, so its results were very similar to the 
aluminum sphere's results off resonance. Waveforms were practically identical for 
the shorter pulse lengths regardless of frequency and composition implying the 
reflections are purely specular. The results clearly demonstrated the dependence of 
pulse length on reflected waveforms. 

D. SELECTION OF STYROFOAM 

Since defining target strength becomes difficult when using elastic spheres, 
which was not the intent of this work, the need existed to do one of two things: (1) 
work at a frequency that gives good replication of the input pulse, or (2) find a 
target that gave a very high impedance mismatch. 



34 




R e I ai 1 ve T i me 



Figure 16. Aluminum Sphere Reflection at Resonance, 82.7 kHz 
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Figure 17. Aluminum Sphere Reflection, 84.7 kHz 
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Figure 18. Aluminum Sphere Reflection, 80.7 kHz 
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Figure 19. Brass Sphere Reflection, 80.7 kHz 
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Figure 20. Brass Sphere Reflection, 82.7 kHz 
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Figure 21. Brass Sphere Reflection, 84.7 kHz 
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Figure 22. Aluminum Sphere Reflection, 82.7 kHz 
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Figure 23. Aluminum Sphere Reflection, 84.7 kHz 

N=2 
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Figure 24. Aluminum Sphere Reflection, 82.7 kHz 
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Figure 25. Aluminum Sphere Reflection, 84.7 kHz 

N=5 
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Figure 26. Aluminum Sphere Reflection, 82.7 kHz 
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Figure 27. Aluminum Sphere Reflection, 84.7 kHz 

N = 15 
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Figure 28. Brass Sphere Reflection, 82.7 kHz 
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Figure 29. Brass Sphere Reflection, 84.7 kHz 
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Figure 30. Brass Sphere Reflection, 82.7 kHz 
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Figure 31. Brass Sphere Reflection, 84.7 kHz 
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Figure 32. Brass Sphere Reflection, 82.7 kHz 
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Figure 33. Brass Sphere Reflection, 84.7 kHz 

N = 15 
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The most important factor in achieving a good reflection is the impedance 
mismatch between the fluid and object. Experimental work done with metal spheres 
by Hamptom [Ref. 5] shows the effect of compostion using several metals. A 
pressure wave in water will see little resistance to a hard rubber sphere because its 
characteristic impedance is near that of the water. An object that gives a high 
impedance mismatch is an air bubble. The impedance ratio of water to air is 
roughly 3000. For the case of an air bubble, the boundary condition between the 
water and air is termed ’’pressure release”, since the pressure transmitted into the air 
is nearly zero [Ref. 2:p. 127]. The difference between a pressure-release sphere and 
a rigid sphere is best seen graphically in Figure 34 [adapted from Ref. 15: p. 359]. 
In the ka region used in the experiments, both the rigid and pressure-release form 
functions will be nearly identical and equal to one. 




Figure 34. Backscattered Form Function of a Rigid and Pressure 

Release Sphere 

A styrofoam ball was used to approximate an air bubble not near resonance. 
The reflected waveforms were generally identical to the input pulse at all 
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frequencies. This made styrofoam an ideal target for this work. Figures 35 
through 40 show reflections from a three-inch diameter styrofoam sphere using 
various pulse lengths at 100 and 106 kHz. 

E. BISTATIC SCATTERING 

A 360 degree pattern of the pressure amplitude scattered from a styrofoam 
sphere was obtained to compare the results with theoretical predictions for a 
perfectly reflecting sphere. Since a sphere is perfectly symmetrical, it was not 
necessary to take the entire 360 degrees of data but only 1 80 degrees of it. 

Theoretical equations have been derived and published by several authors that 
predict scattering of acoustic plane waves from a sphere [Ref. 16,17]. They tend to 
be quite complicated functions involving sums of Hankel functions. Figure 41 was 
adapted from Medwin and Clay [Ref. 16:pp. 192-193] and shows computed 
scattering patterns in relative amplitude for a perfectly reflecting sphere in water 
for several ka values. The incident plane wave is arriving from the right. The 
patterns show that for increasing ka the major forward lobe gets larger in amplitude 
and narrower in width. More structure, i.e., secondary lobes, develop as ka 
increases. 

Bistatic data were obtained for a ka =10 using a three-inch diameter styrofoam 
sphere and are shown in Figure 42. The incident wave arrives from the bottom. 
Frequency was 61.5 kHz and an eight cycle pulse was used. The graphed data 
revealed some structure in the forward scattering region and a large forward 
scattering lobe. Comparison with the theoretical prediction for a rigid sphere 
(dashed line in Figure 42) show that the salient features of the scattered pattern 
were duplicated. 
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Figure 35. Styrofoam Sphere Reflection, 100 kHz 
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Figure 36. Styrofoam Sphere Reflection, 106 kHz 

N=5 
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Figure 37. Styrofoam Sphere Reflection, 100 kHz 

N = 10 




Figure 38. Styrofoam Sphere Reflection, 106 kHz 

N = 10 
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Figure 39. Stvrofoam Sphere Reflection, 100 kHz 

N = 20 




Figure 40. Styrofoam Sphere Reflection, 106 kHz 
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Figure 41. Theoretical Scattering Patterns for 

Rigid Spheres in Water 
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VI. RTSTATTC SCATTERING OF A STYROFOAM CYLINDER 



A styrofoam cylinder of diameter 2.5 cm and length 20 cm with flat ends was 
used for bistatic scattering measurements. The particular diameter was chosen 
because the styrofoam was commercially available in that size. The length was 
selected using two criteria: (1) the length was small enough so that the receiver 
could be placed in the far field of the sound reflected from the cylinder and (2) it 
was large enough to clearly present a cylindrical shape. 

Since the F41 was not omnidirectional, there was concern of approximating a 
plane wave at the target. The lower the frequency, the better the pressure 
distribution across the target because of the larger major lobe, but the wavelength- 
to-target length ratio would be less and the signal-to-noise ratio also became poorer 
at lower frequencies. A frequency of 70 kHz was chosen at a source-to-target 
distance of 125 cm. Figure 43 shows the incident pressure distribution as a function 
of distance off the axis at a distance of 125 cm. The lower curve is for the 
monostatic case where the receiver shadowed some of the incident energy with a 
receiver-to-target center distance of 50 cm. A two percent pressure amplitude 
difference existed for the shadowed case. An eight percent decrease in pressure 
amplitude existed when not shadowed. 

The cylinder was mounted to a 1/8-inch diameter metal tube which had a 2 X 
1/2 X 1/1 6-inch flat plate epoxied to it to provide a means of rotating the cylinder. 
The flat plate was forced into the interior of the cylinder and provided a snug fit. 
The effect of the plate on the reflections was assumed to be neglible since it was 
completely surrounded by the styrofoam. The tube was as small as possible to 
minimize reflections. A string with weights was attached to keep the st>Tofoam 



51 




Distance from Center 

Figure 43. Pressure Distribution at 70 Khz 
cylinder submerged. The thin tube was soldered into a 3/8-nch diameter tube 
which supported the entire arrangement. The cylinder was level before entering the 
water and appeared to remain that way when in the water. Figure 44 illustrates the 
support setup . 

Data were recorded for nine bistatic angles. Time discrimination of the 
reflected pulse was used whenever possible. The source-to-receiver distance was 
125 + 2 cm in all cases with the F41 being driven at 177 V rms. The Ithaco 
preamplifier setting was 20X with a 3-300 kHz bandwidth. Receiver-to-cylinder 
center distance varied, but the receiver was always in the far field of the cylinder. 
All data were extrapolated to one meter from the target center for comparison. A 
20 cycle pulse was used. 

The cylinder was initially aligned broadside to the source by eye to within an 
uncertainty of about + 4 degrees. Assuming symmetry, data were taken for 180 
degrees of rotation and several points were taken beyond t hat limit. Figure 45 
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Figure 44. Support Platform Setup 



illustrates the angles measured in the experiments. The varied angle. A, is the angle 
made between the cylinder's normal and the line connecting the center of the target 
to the source axis. Angle B is the bistatic angle. 



Source 




Figure 45. Measured Angles in Experiment 
The results of the experiment are shown in Figures 46 to 55. In most cases, 
strong peaks were observed when angle A was one-half the bistatic angle and 
secondary peaks with smaller amplitude occurred 90 degrees from that. This 
geometry seems to indicate specular reflection off the long cylinder's surface and its 
flat ends respectively. As the bistatic angle increases the secondary peak broadens 
and eventually disappears. Well into the forward scattering region, greater than 
150 degrees of bistatic angle, the data became less predictable The distinction of the 
primary and secondary peaks became less clear with additional peaks showing up. 
At a bistatic angle of 180 degrees no distinctive peaks were present and amplitude 
became very large in comparison to other bistatic angles. 
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Figure 46. Styrofoam Cylinder, Backscattered, 70 kHz 
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Figure 47. Styrofoam Cylinder, Bistatic Angle 10 Degrees, 70 kHz 
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Figure 48. Styrofoam Cylinder, Bistatic Angle 20 Degrees, 70 kHz 
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Figure 49. Styrofoam Cylinder, Bistatic Angle 30 Degrees, 70 kHz 
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Figure 50. Styrofoam Cylinder, Bistatic Angle 50 Degrees, 70 kHz 
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Figure 51. Styrofoam Cylinder, Bistatic Angle 90 Degrees, 70 kHz 
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Figure 52. Styrofoam Cylinder, Bistatic Angle 120 Degrees, 70 kHz 
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Figure 53. Styrofoam Cylinder, Bistatic Angle 150 Degrees, 70 kHz 
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Figure 54. Styrofoam Cylinder, Bistatic Angle 165 Degrees, 70 kHz 
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Figure 55. Styrofoam Cylinder, Bistatic Angle 180 Degrees, 70 kHz 



VII. CONCLUSION AND RECOMMENDATIONS 



A data acquisition system using the Tektronix 2430 digital oscilloscope was 
implemented to extract scattered waveforms in a noise-filled water environment. 
The system was verified to work and reflected waveforms were collected. 
Preliminary work revealed that the relative movement of the source and receiver 
between collection of waveforms was critical in obtaining accurate data using a 
digital sampling system with a subtraction routine. Data collection was slow and 
tedious. 

Several aspects influencing scattered waveforms, including target composition, 
frequency, and pulse length were observed and recorded. The effects of these on 
the waveforms of the reflected pulse produce difficutlies in measuring the intensity 
of pulses with waveforms not identical of the incident pulse. 

Styrofoam targets were used to circumvent the problem of waveform 
dependence on frequency, pulse length, and composition. Bistatic measurements 
were taken using a sphere and a cylinder as targets. The data from the sphere 
compared well with theoretical predictions for a perfectly reflecting sphere. 
Readily available data were not found for comparison with the results obtained for 
the cylinder. The goal of matching wavelength-to-target length ratios observed in 
oceanic conditions was not done in this work. The USRD type E27 transducer 
would be an excellent source for this application since it has a frequency range 100 
to 500 kHz. Its small size makes its beamwidths very large which is beneficial when 
receiver-target distance is limited by tank dimensions. The E27 is also very 
lightweight and easily moved. The effect of the jitter would have to be evaluated at 
the higher frequencies. 
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The bistatic radar cross-section theorem could not be compared to the data 
collected in this report. The wavelengths used in the experiments are not close to 
the "vanishing limit" as stated in the theorem. One common trend between the 
bistatic data and the theorem was that the data became less predictable as the bistatic 
angle approached 180 degrees. 

For future work, a physical structure must be made to facilitate positioning of 
the source, target and receiver in three dimensions so that data collection can be 
completely automated and accuracy improved. The rotation of supports should 
have a digital device to allow control by the computer. The design and construction 
of such a system would not be a trivial project and would require a great amount of 
ingenuity and experience with electronic devices. The payoff of having such a 
system would be great to students studying underwater acoustics and in particular to 
anyone studying scattering effects. Once a positioning system is built and 
completely automated, data collection will become simpler and more time can be 
devoted to analyzing it. 
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APPENDIX A 

COMPUTER LISTING OF "TEX" 



PROGRAM TEX 



DATA AQUISITION FOR BI-STATIC TARGET STRENGTHS 
USING THE TEKTRONIX Z430 DIGITAL OSCILLOSCOPE 

THIS PROGRAM TRANSFERS WAVEFORM DATA FROM THE TEX Z430 
TO THE COMPUTER WHERE IT IS CONVERTED TO ACTUAL VOLTAGE 
VALUES AND STORED ON A FLEXIBLE DISK. INPUT FROM THE 
USER REQUIRES FILE NAME, DESCRIPTION. AND CONTINUATION 
RESPONSES. PREVIOUS SETUP OF THE TEX Z430 IS REQUIRED. 

A THOROUGH KNOWLEDGE OF THE GPIB CODES IS NECESSARY. SEE 
THE OPERATOR’S MANUAL AND THE INTERFACING MANUAL FOR 
THE TEX Z430. 



INITIALIZATION 

A J, A ^ A.A. Xa A A. A.A A A A A A A AA *1 A » A A fc* At Al, A A J ^ J J. A A. A > >1 -tA M k. A A. A± A k -lA A* Aa Aa M. A A M. A ±A A A. MM AA A-M. AA MA AA JA AA AA Al ^4 At AA lA 

TtitTrit t rTtiATr TT n ii it i! z rTt ti iiiTirT rT r i T i T tri iu tr it 

tT TT T ? T f • I TT ▼T TT TT TIT? M M IfTt TT TT TT 'TT ^T TT tT TT TT TT TT TT TT Tf T V TT TT TT TTTfTTTTTTif II M If fl fl if 1 1 TTTTTTTTTT Tr TT TT ?f Tf TT 9 f PI f I Tf f f Tf 

PRINTER IS CRT 
OPTION BASE 0 

DIM Vdata$[ 151 ,W$140] INAME OF DATA FILE AND DESCRIPTION 
DIM A$[Z00] ,B$t 5000] ,D$I 5] ! PRE-WAVEFORM AND WAVEFORM STRINGS 
DIM V( 10Z3),Tim( 10Z3),VoU( 10Z3)!DATA ARRAYS 
DIM X$[ 10] ,Ys$t 10] ,Vs$I 10] ,Ts$[ 10] 

DIM An$[ 1 ] 



* 4 A A A 41 A A. A A A A A A A A A A Ml » h AA » M MM AM AAAA AA AAAA A i MM Ml MA Ml A A A A A A AA AAA* k M A A * ■ * 1 Ml Ml A A MM A A li A I AM II AA M M A A AA At MM AA AA A A AA AA A A AA AA 

ffTTirilir n n iiTTTrTrTtTtTtTritlTirTfirTTTTTITTiTTrTTiTTiItltir n nTTlTItltirTTn ITirilTrTTTT 

TTtT TT T ? TT TT TT T? 1 r “ f " T TT TTTT ^T ^TTT TT TT^T TTTTTT TT TTTT TT TT TTTT TTTT TT^T TT TT TT TTTTTTTTTT TTTT TT TT TT TT^T^TTtT TT T ? TT TT TT TT TT^^TTT 

TEKTRONIX SETUP 



OUTPUT 704; "PATH OFF" ! SHORTENS THE DATA STRING SENT BY THE Z430 



INPUT FROM USER 

AA AA AA AA Al AA AA AA AA A> AA A A A A Al AA AA »a ai aa AJ MM MM JA A* AA Mm AA AA AA A A A A AA AJ AA AA AA AA AA AA AA AA AA A ■ A A AA A A AA AA AA AA AA AA AA A> AA AA A A A A A A A. A A.A A-A- AA -A^ 

ITlTirTI LI Tiltli it u tTlTgiT “ ^i rnuttiiti ITui r^^^P II I ! nTT itTTItItT r iTTi iT » 

TTTT^X ?f f I TTTT TT TT M ft TT TT TT TT fT ff ff TF FPTTTT TTTT TT f T ff TTTT TT TTTT TT f f f f ? f TTTTTT TTTTTT TT^T TT TTTTTT TTTTTT TTTT TT^T TT TT TT TTTT TTTTtTTT 



INPUT "NAME OF THE DATA FILE?" . VdataS 
INPUT "DOES FILE ALREADY EXI ST?( Y=YES ) " , An$ 

IF An$="Y" THEN PURGE VdataS 

INPUT "INPUT FILE DESCRIPTION?( 40 CHARACTERS MAX)",W$ 

BEEP 

PRINT "ALLOWING DATA DISK INSERTION. PRESS CONTINUE TO RESUME." 
PAUSE 
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500 

510 

5Z0 

530 

540 

550 

5B0 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 



CREftTE BDftT VdataS , 1 , 1 7000 
ASSIGN §Dfile TO VdataS 
PRINT 



AM. » A Ah. A*. A* ->.4. AJ. AM M K ^ A MA ^ A At A.A. A A AM M M A A MM A% AA MA MA A. M. MA A. A A A A A A. A A A A. A Aa A^ AA A. a A.A AM aA Am. a a a. A At. >> A A am a A A a A A A A Aa MA mA A.A AA AM A A A A 

ittiTT tTTlttTtlt it till tttru 

TV TT W Tt TT W W W TT W W WTT 

OBTAIN WAVEFORM 



PRINT "OBTAIN DESIRED WAVEFORM ONTO TEKTRONIX" 
PRINT "PRESS CONTINUE TO RESUME" 

PAUSE 



A A A A XX XA A A XX AX A_4 ^X ^X^XXA XX^X^^,^XXX^X^X ^X AJ^AA^^^AA^X AA 

TT IT TrTT TT Tl II fl Tl fl II II II fl U ti 11 TT xi ri ti ii it it n itTTTT^I fl II TTTT T TTi If M II TT Tt if It 11 ii i_i II Tl II II ^1 T Tri TT li iT Ti TT Ti TTii II fi t l II 

W WW WaWt T r T r w w xw W TT 

RECEIVING PRE-WAVEFORM DATA AND THE WAVEFORM 



OUTPUT 704;"WFM? 
ENTER 704 ;X$ 
OUTPUT 704; "WFM? 
ENTER 704;Ys$ 
OUTPUT 704; "CHI 7 
ENTER 704;Vs$ 
OUTPUT 704; "HOR? 
ENTER 704;Ts$ 

! 



XIN"!TIME INCREMENT BETWEEN POINTS 
YMU"! VOLTAGE DIGITIZING INCREMENT 
VOL" IVOLTS/DIV 
BSE"!B TIME SCALE 



I 



OUTPUT 704;"CURV?" 
ENTER 704 ;B$ 



AA AM M A MM A k *1 MM a ^ 41 *> Mi A A AM A.1 A A A A A> >■ AA **AA*A M A A A MM AM A A AM At. At A>, 1. M AM MM A A 4i >A AA AA AA m a MM Mm m M A a A A A A ^ A AM a A AA-AA AA. AA AA AA AA AA AA AM A A M k. 

tttTttiitr^irtrH it ttl^ttutt ttuiZ ni t uiititt t tuuuuTt 

TT^^rTWTTTTTTTTT^TTTTFTTTVTT^T^T^T W^WTTtTTTWT^TTTTWT^TTWTTTTT TVW^T TV W^VW^V^VW^V^V^TT^T^T ^TTT ^WT^WTTTWTWTTWx TT^T 

SEPARATION OF ASCII STRING INTO DATA POINTS 

A A AA AA A A AAAAAA AA AA AA AA AA AAAA AA A A AA AA A> AA A> MM A A A A AA AM A A AM AA A^ fcA 4A AA tA AA A k A A AM M A AA AA A.A A * M A XX MM A A A A AX A< M A m A MM At AM ll AA At A A A-A A A ^A 

TV I ? wf ^TTT f I wf ^V^V M fl II II ff If ff fl 11 ft f f ft TV Tf If if ^TTT^T^T If T f f f f f TTTVTT^V^VTvTVTT Tf f f If ^V^^VTTV ¥f If II II II y f 3 f If ff TT WTV 



N=LEN( B$) 

D$="" 

I nc=0 
! 

FOR 1=1 TO N 
IF B$n ,n=" THEN 
V( Inc)=VRL(D$> 

Inc=I nc+ 1 
D$="" 

GOTO 950 
END IF 

D$=D$&B$[ I , I ] 

NEXT I 



CONVERSION TO REAL-TIME DATA 

«« «.! AA Jl. m A .A X Xt A.M. t. A 1.1. »« i ± «t «« »« «» M A A ^ A A AA AA A.A A1 .lA AA AA A. A. AA At «A At ■■ AA AA AA A« AA AAAA »A A« AA AA A. .a AA AA AA 

itfl n TfiT li Tl ITir it iT ii ITTT 11 ^ ** tlTl i t^*» ** TI tTTT TTTttT tTTt ttTTtT TT T* tT TTTT TT TT TrtTTrn n TTI 1 ] jl Tiff nn It Tr ^TrTTTT TT TT 

TTTT TrTT If fl M iV II lift If f ? ^V^TTT If Ff ff if If fl TT^TTT^T^TTT i I • I TV TV Tf 1 ? TV ^T TrTT f f Tf TT^T^TTT TfTTTfTI If fl II II ff If fl TT^T^T^T f r f I f 1 
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1010 ' X$=TIME INCREMENT BETWEEN POINTS DIV/50 

1020 ! YS$=VOLTS BETWEEN POINTS DIV/Z5 

1030 I VS$=VOLTS/DIV 

1040 I TS$=B TRIGGER TIME SCALE 

1050 I 

1060 DeItin=VAL(X$) 

1070 Delvol=VAL( Ys$) 

1080 Vscal e=VAL( Vs$ > 

1090 Tscale=VAL( Ts$ ) 

1 100 ! 

1110 OUTPUT @Df ile;W$.Tscale,Vscale 
1120 FOR 1=0 TO 1023 
1130 TiPi( I ) = I*Deltim 
1140 VoI-t( I > = V( I )*Delvol 
1 150 NEXT I 
1160 ! 

1170 PRINT "WRITING TO FILE" 

1180 OUTPUT @Dfile;Tin(*),VoH(*) 

1 190 ASSIGN ©Df ile TO * 

1200 BEEP 
1210 > 

1220 I tiltltltttrttttttt ttttttltttttltttttttltttlttt'tfttttttttttlttttlTt'tttTtTtttl tttltttttttttttttttlitttltltttllTITtttlt'ttlt II II tl tl 

I ^ ^ T f f T TT II f 1 T^TTTT TT TTTT if IP M rw tT f| II f| If If tt f I ?F TT Fffl II wf ff “T T? f T Tf “ f “ f TT “ f ^ TT'^T^ TT tT 

1230 ! 

1240 INPUT "DO YOU WANT ANOTHER WAVEFORM^ Y=YES",An$ 

1250 IF An$="Y" THEN GOTO 430 
1260 ! 

1270 PRINT "PROGRAM TERMINATED" 

1280 END 
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Z0 
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1 10 

1 20 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

2 30 

240 

250 

260 
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300 
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370 
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390 
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420 

430 

440 

450 



APPENDIX B 

COMPUTER LISTING OF "PLOT" 

PROGRAM: PLOT 

THIS PROGRAM TAKES TUO FORMATTED BOAT FILES AND SUBTRACTS 
ONE FROM THE OTHER. IT WAS DESIGNED FOR USE USING WAVEFORMS 
OBTAINED ON THE TEKTRONIX 2430 DIGITAL OSCOPE. DATA IS TAKEN 
FROM THE FILES IN THE FOLLOWING ORDER. 

W$ = FILE DESCRIPTION STRING 

TSCALE = TIME SCALE OF THE SCOPE ( SEC/DI V) 

VSCALE = VOLTAGE SCALE OF THE SCOPE (VOLT/DIV) 
T(0-1023) = TIME ARRAY (SEC) 

V(0-1023) = VOLTS ARRAY (VOLTS) 

THE SUBTRACTION RESULTS ARE THEN ABLED TO BE PLOTTED AS HARD 
COPY ON THE THINK JET OR PEN PLOTTER. A STANDARD OPTION IS 
INCORPORATED FOR EASE WHEN USING THE PROGRAM FREQUENTLY WITH 
THE SAME PARAMETERS. THE PROGRAM ALLOWS FLEXIBILITY IN 
CHANGING THE PLOT’S PARAMTERS BEFORE THE HARD COPY. 

AM -a ■ AM M. m. Aa AM AA fc A AM M t M. M ai Al. t A M. A AM A.M A.A AAAAAa m a A A MM MM AAAAkM AM AM A A A A A A A A A A A A A A A A M A M A A A AA M A AM AM 4.A M A. AA 

IT Tl Tt II n tl TrTTTtri^TT TTTTTtTT^TTt T ^TT^ t7 TTTT TT TT TT H t f Tt tt IT TT TT tT tT Tttt^TTTT T T till 11 IT 11 21 II IT l Tllirir II Ti 11 H II TTiTTr 

f 1 TT If TT f1 iJTf ▼T f f f ? Tf ‘TT Tf f f f f T f TTTT‘ f f 1 f f T ?f TTrT Tr ff TT ▼T TT TT TT ffl PI fl fl II II ff TT II •! M fl M TfTF ff IF ff ▼TTT 

INITIALIZATION 



OPTION BASE 0 

DI M X_da't a( 1023), Y_da'ta( 1023), Y_ax is_uni'ts$[ 32 1 
DIM Head_lbl$[ 60] , Y_ax is_lbl $[ 60] , X_ax i s_l bl $l 32 ] 

DIM X_ax is_uni t s$[ 32 ] , Cl ear_crt $[ 2 ] , C$[ 2 ] , Sub_head_l bl $[ 32 ] 
DIM X_char$[8] ,Vdata$[ 15] ,W$[40] .VdataZSl 15] ,Y_char$[8] 

DIM XI ( 1023) ,Y1( 1023) ,X2( 1023) ,Y2( 1023) 
Clear_crt$=CHR$(255)&CHR$(75) 

C$=CHR$( Z55)&''K" 




GOSUB Load_file 

INPUT "DO YOU WISH TO USE THE STANDARD PARAMETERS? Y=YES" , Dec2$ 
IF Dec2$="Y" THEN 



GOSUB Stand_parans 
GOTO 420 
END IF 

GOSUB Init ial_parar»is 
A$="N" 

GOSUB Plot_data 

INPUT "WISH TO CHANGE CURRENT GRAPHICS PARAMETERS? Y=YES",Dec2$ 
IF Dec2$="Y" THEN GOSUB Change_params 
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460 INPUT "ANOTHER DATA FILE PLOTTED ON SAME HARD COPY GRAPH?" , Decb$ 
470 UHILE DecS$="Y" 

480 GOSUB Load_file 

490 INPUT "LINE TYPE! 1 =SOLI D. 4=DASH . 5=BIG DASH . G=DASH-DOT" , L i ne_num 

500 INPUT "PEN NUMBER 1 TO 8" ,Pen_num 

510 GOSUB Sane_graph 

5Z0 GOTO 460 

530 END UHILE 

540 INPUT "DO YOU WISH TO CHANGE ANY MORE PARAMTERS? Y=YES".DecZ$ 

550 IF DecZ$="Y" THEN GOTO 450 

560 INPUT "DO YOU WISH TO DUMP GRAPHICS? Y=YES",A$ 

570 IF A$»"Y" THEN GOSUB Plot.data 

580 INPUT "DO YOU WANT TO PLOT A DIFFERENT GRAPH? Y=YES",Dec7$ 

590 IF Dec7S="Y" THEN GOTO 350 

600 GOSUB End_progran 

610 ! 

630 ! STANDARD PARAMETERS 

650 ! 

660 Siand_params : • 

670 Head_lbl$=U$ 

680 Sun_head_l bl $=Vdat a$ 

690 X_char$*"K" 

700 Y_char$="DD.DD" 

710 Y_ax 1 5_1 bl $=" Vol t s" 

720 X_ax i5_l bl J *" Rel at 1 ve Tine" 

730 Y_ax i5_unit5$="Vol ts" 

740 X_ax is_units$="nicroseconds" 

750 Left=0 

760 Right= 1 9»Tscal e 

770 Botton*=-Vscale*4 

780 Top*Vscale*4 

790 X_t ic_sp*Tscale/Z 

800 Y_t 1 c_sp=Vscal e/5 

810 X_naj_ct=Z 

8Z0 Y_naj_ct=5 

830 Ma]_tic*=Z 

840 Xgridt ic*Tscale 

850 Ygridt ic»Vscal e 

860 Xgrldnaj=*l 

870 Ygridna]“l 

880 Majgrxd_5Z*Z 

890 Y_incnt=Z * Vscal e 

900 X_incnt*4*Tscale 

910 Line_nun*=l 

9Z0 RETURN 

930 ! 

M I ** k A ^ * M.K M.M. A ^ AA. AA. Aa a AA a* A* «A AA a-A. A A A A A-». AA aA AAAA A- a AA *■ ^ A ^ AA. aA A A AA. A A. A A- A^A^-4-4- A ». AA. 

940 ! ##ggggggggggtig#gggggg#lMtlUtttll tf 4g4ggttw 4 U4 1t tttttt4tgtt4i4U»t n l^ 

950 ! INITIAL PARAMETERS 
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970 
980 
990 
1000 
1010 
10Z0 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1 100 
1 1 10 
1 1Z0 
1 130 
1 140 
1 150 
1160 
1 170 
1 171 
1 180 
1 190 
1 Z00 
1Z10 
1ZZ0 
1Z30 
1Z40 
1Z50 
1ZB0 
1Z70 
1Z80 
1Z90 
1 300 
1301 
1 310 
1 3Z0 
1330 
1 340 
1350 
1360 
1370 
1 380 
1390 
1400 
1410 
14Z0 



Initial _parans: ! 

X_t 1 c_so=0 
Y_t 1 c_sp=0 
X_maj_ct= 1 
Y_maj_ct= 1 
f1aj_t ic=Z 
Xgr idt ic=0 
Ygridt ic=0 
Xgr idnaj = 1 
Ygr idna] = 1 
Maj grid_sz = Z 

INPUT "GRAPH HEADING LABEL . 3Z CHAR. hAX" ,Head_lbl$ 

INPUT "GRAPH SUB LABEL. 3Z CHAR. MAX" , Sub_head_lbl $ 

INPUT "IMAGE FOR X AXIS NUMERICS, 8 CHAR. MAX",X_char$ 
INPUT "IMAGE FOR Y AXIS NUMERICS, 8 CHAR. MAX",Y_char$ 
INPUT "Y AXIS LhBEL,3Z CHAR. MAX" . Y_ax i s_lbl $ 

INPUT "X AXIS LABEL, 3Z CHAR. MAX" , X_ax i s_l bl $ 

INPUT "Y AXIS UNITS, 3Z CHAR. MAX" , Y_ax is_units$ 

INPUT "X AXIS UNITS, 3Z CHAR. MAX" , X_ax i s_units$ 

DISP "UINDOU LOOKED AT .LEFT . RIGHT . BOTTOM , TOP" ; 

INPUT Lef t , Right , Bot ton , Top 

INPUT "X AXIS TICK SPACI NG , DEF = 0( NO TI CKS ) " . X_t ic_sp 
INPUT "Y AXIS TICK SPAC I NG . DEF=0( NO TICKS >" ,Y_tic_sp 
INPUT "NUM TIC INTVLS BTUN MAJ TICS X-DIR,DEF=1 " ,X_naj_ct 
INPUT "NUM TIC INTVLS BTUN MAJ TICS Y-DI R , DEF» 1" , Y_naj _ct 
INPUT "SIZE OF MAJOR TICKS. DEF=Z" ,Maj_tic 
INPUT "GRID X-TIC SPACING ,DEF=0< NO TICKS)" .Xgridtic 
INPUT "GRID Y-TIC SPACI NG . DEF=0( NO TICK S )". Ygridt ic 
INPUT "NUM TICKS BTUN MAJOR TICKS Y-DI R , DEF= 1" , Ygr idnaj 
INPUT "NUM TICKS BTUN MAJOR TICKS X-DI R . DEF= 1 ", Xgr idnaj 
INPUT "SIZE OF MAJOR GRID TI CKS . DEF=Z" . Maj grid_sz 
INPUT "Y-AXIS INCREMENT" ,Y_incnt 
INPUT "X-AXIS INCREMENT" ,X_incnt 

DISP "LINE TYPE( 1=S0LID.4=DASH,5=BIG DASH.6=DASH-D0T>" ; 
INPUT Line_nun 
! 

RETURN 

! 




I 

Load_file: ! 

! 

INPUT "INPUT NAME OF DATA FILE 1",Vdata$ 
INPUT "INPUT NAME OF DATA FILE Z ".VdataZS 
ASSIGN @Df1 TO VdataS 
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1430 ENTER @Df 1 ; U$ , Tscal e . Vscal e . X 1 ( » ) , Y 1 ( • ) 

1440 ASSIGN eOfl TO • 

1450 ASSIGN @DfZ TO VdataZ$ 

1460 ENTER @Df Z ; U$ , Tscale . Vscale , XZ ( * ) . YZ( ♦ ) 

1470 ASSIGN @DfZ TO * 

1480 MAT X_daia= XI 

1490 MAT Y_daia= Yl-YZ 

1500 PRINT "DATA FILES READ" 

1510 BEEP 
15Z0 RETURN 
1530 ! 

H ii vj I XT TTTxTTn p II H " w n iitm IiTiit n it it tti t ^ ^^^””""""*4n irnHiiffHn7Tiiiinuiin httii n irn 

1550 ! PLOT DATA 

t r* /^ a A A A A* AA AA A-A a a a a AA A A A ^ A A A A A A Aa A A A A a A A A AA Aa M A JA Aa AA A* AA AA Aa Aa A* *A AJ ^A Aa A* A^ AA AA AJ AAAA^ A A A-A A-A A.A A A AA A_A A^ A^ A A AA A A A A A A A A 

SriW 1 1 I H TTir u nTT II tiii4r« II H iiiiiiairirTr cirnnnitirirnir^ g p n fl ircn i rii H T TT TiTTrTi n ir n tt ii n n Tiiriiiriiiiir 
1570 ! 

1580 Plo‘k_da'ta: ! 

1590 ! 

1600 IF Fl$="Y" THEN 

1610 OUTPUT KBD;Clear_crt$; 

16Z0 GINIT 

1630 PLOTTER IS CRT , "INTERNRL" 

1640 GOTO 1770 

1650 END IF 

1660 INPUT "IS THIS TO BE A HARDCOPY PLOT? Y»YES",Decl$ 

1670 IF Decl$="Y" THEN 
1680 GINIT 

1690 PLOTTER IS 705,"HPGL" 

1700 GOTO 1770 

1710 ELSE 

17Z0 OUTPUT KBD;Clear_crt$; 

1730 GINIT 

1740 PLOTTER IS CRT ," INTERNAL" 

1750 GOTO 1770 

1760 END IF 
1770 GRAPHICS ON 
1780 PEN 1 

1790 X_gdu_nax=100*MAX( 1 .RATIO) 

1800 Y_gdu_nax-100«MAX( 1 , 1/RATIO) 

1810 LORG 5 
1820 CSI2E 4 

1830 MOVE . 48»X odu wax..93»Y odu wax 
1840 LABEL Head_lbl$ 

1850 MOVE . 48»X odu wax..89*Y odu nax 

I860 CSIZE 3 

1870 LABEL Sub_head_lbl$ 

1880 DEG 
1890 LOIR 90 
1900 CSIZE 4 
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1910 
19Z0 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
Z000 
Z010 
Z0Z0 
Z030 
Z040 
Z0S0 
Z060 
2070 
Z080 
2090 
2100 
Z1 10 
Z 120 
ZI30 
Z 140 
Z150 
Z160 
Z170 
Z180 
Z 190 
ZZ00 
Z210 
ZZZ0 
ZZ30 
ZZ40 
ZZ50 
ZZB0 
ZZ70 
ZZ80 
ZZ90 
Z300 
Z310 
Z3Z0 
Z330 
Z340 



MOVE .04*X gdu nax..55*Y odu nax 
LABEL Y_axis_lbl$ 

LORG 5 
LDIR 0 

MOVE . 48*X_gdu_max , . 1 5*Y_gdu_nax 
CSIZE 4 

LABEL X_axis_lbl$ 

MOVE . 48*X_gdu_nax , . l9*Y_gdu_nax 
CSIZE 3 

LABEL •■(";X_axis_units$;")" 

MOVE .70»X gdu wax..74»Y odu wax 

VIEUPORT . 1 5*X_gdu_nax , . 85*X_gdu_max , . Z5* Y_gdu_max , . 8G*Y_gdu_nax 
UINDOU Lef t , Right , Bottom , Top 

AXES X_t ic_sp , Y_t ic_sp , Lef t .Bottom, X_maj _ct , Y_maj_ct , Maj _t i c 
AXES X_tic_sp,Y_tic_ 5 p, Right ,Top,X_maj_ct ,Y_maj_ct ,Maj_tic 
• GRID Xgr idt ic , Ygndt ic .Left .Bottom , Xgr idmaj , Ygndmaj .Majgr id_sz 
CLIP OFF 
CSIZE 3. .75 
LORG G 

FOR I=Left TO Right STEP X_incmt ! LABELI NG OF X AXIS 
MOVE I , ( Bot t om- . 04» ( Top ) ) 

LABEL USING X_char$ ; I / I . E-G 
NEXT I 
LORG 8 

FOR I=Bottom TO Top STEP Y_incmt ! LABELING OF Y AXIS 
MOVE ( Lef t-.4*Tscale) ,I 
LABEL USING Y_char$;I 
NEXT I 
PENUP 
PEN 2 

LINE TYPE Line_num 
CLIP ON 

FOR 1=0 TO 10231PLOTTING OF DATA 
PLOT X_data( I ) ,Y_data( I ) 

NEXT I 
PENUP 
CLIP OFF 

IF A$="Y" THEN DUMP GRAPHICS «701 
RETURN 




Z 350 Same_graph: ! 

Z3G0 PEN Pen_num 
Z 370 LINE TYPE Line_num 
Z380 CLIP ON 
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Z390 FOR 1=0 ro 1023 

Z400 PLOT X_data( I ) , Y_data( I ) 

Z4I0 NEXT I 
Z4Z0 CLIP OFF 
Z430 PENUP 
Z440 RETURN 
Z450 ! 

Z4b0 I gw4M4444T444H4g 44w 4tt44W M t HH w4HHUtt4tg 41t > n »44MTff4ttttC 4 4Utt4Mt41ffl4#4»» 4 ftH 
Z470 ! CHANGE GRAPHICS PARAMETERS 

^ I AA 1* *.i AA A k. AA AA AA Aa Aa A A AA A A A A AA ^A AAAA -*-A *-A »A «A A. * At ■■* ** *> AB. A ^ M. A JA *.i AA A. A Aa A A A ■< A A AA A A A.A AA AA .AA. A-^ AA. A-A ^AA AA ^AA. AA- A-A A-A A_A_ A-4^ A-A 

/awi/l I 11 H Ft"TyTT iT TF n TT P II H 11 iT IT Tl ii rl o rf Ti IT it Ti ii n H U FiiT Ti ii ft Tt llnll lt UrT xI tlH tl lllltlHTYiinYTriWTIHIrTTHn TrTT TTTT 

Z490 • 



Z500 Change_parafns : • 



Z5I0 
Z5Z0 
Z530 
Z540 
Z550 
Z5G0 
Z570 
2580 
Z590 
ZG00 
ZB 10 
ZGZ0 
ZG30 
ZG40 
ZF:50 
2GG0 
Z670 
ZG80 
ZB90 
Z700 
Z710 
Z7Z0 
Z730 
Z740 
Z750 
2760 
2770 
Z780 
2790 
2800 
2810 
Z8Z0 
Z830 
Z840 
Z850 
Z8G0 



OUTPUT 
PRINT ■ 
PRINT ■ 
PRINT ■ 
PRINT ■ 
PRINT • 
PRINT ■ 
PRINT ■ 
PRINT ■ 
PRINT • 
PRINT • 
PRINT ■ 
PRINT • 
PRINT • 
PRINT ■ 
PRINT ■ 
P<=INT ■ 
PRINT ■ 
PRINT ■ 
PRINT ■ 
PRINT ■ 
PRINT ■ 
PRINT ■ 
PRINT " 
I 



KBD; Clear_crt$; 

SELECT FROM THE FOLLOWING TO CHANGE" 
GRAPH HEADING" 

Y AXIS LABEL" 

X AXIS 

Y AXIS 
X AXIS 
WINDOW 
X AXIS 

Y AXIS 



CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 

CASE 



1 

Z 

3 

4 

5 
B 

7 

8 
9 
10 ; 
1 1 ; 
1Z: 
1 3: 
14; 
15; 
IB: 
17; 
18: 
19; 
20 : 
Z1 : 
ZZ: 



LABEL" 

UNITS" 

UNITS" 

UNITS" 

TICK SPACING” 

TICK SPACING" 

TICK INTVL BTWN MAJOR 
TICK INTVL BTWN MAJOR 
SI7E OF MAJOR TICKS" 

X TICK GRID SPACING" 

Y TICK GRID SPACING" 
GRID TICK INTVL BTWN 
GRID TICK INTVL BTWN 
MAJOR GRID TICK SIZE" 

Y INCREMENT" 

X INCREMENT" 

LINE TYPE" 

GRAPH HEADING SUB LABEL" 
X AXIS NUMERIC IMAGE" 

Y AXIS NUMERIC IMAGE" 



TICKS IN X DIRECTION- 
TICKS IN Y DIRECTION- 



MAJOR 

MAJOR 



TICKS IN X DIRECTION" 
TICKS IN Y DIRECTION" 



INPUT "DO YOU DESIRE TO CHANGE ANY PARAMETERS? Y-YES" ,Dec4$ 
IF Dec4$="Y" THEN 

INPUT "ENTER PARAMETER NUMBER" , Change_nun 
ELSE 

GOTO 35B0 
END IF 

SELECT Change_nun 
CASE 1 

PRINT "CURRENT HEADING LABEL" ,Head_lbl$ 

INPUT "NEW HEADING LABEL" ,Head_lbl$ 

CASE Z 
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Z870 

Z880 

Z830 

Z900 

Z910 

Z9Z0 

Z930 

Z940 

Z950 

Z9G0 

Z970 

Z980 

Z990 

3<m> 

3010 

30Z0 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

31 10 

31Z0 

3130 

3140 

3150 

3160 

3170 

3180 

3130 

3Z00 

3Z10 

3ZZ0 

3Z30 

3240 

3250 

3260 

3270 

3280 

3Z90 

3300 

3310 

3320 

3330 

3340 



PRINT "CURRENT Y AXIS LABEL" . Y_ax i s_l bl $ 

INPUT "NEU Y AXIS LABEL" , Y_ax i s_l bl $ 

CASE 3 

PRINT "CURRENT X AXIS LABEL" . X_ax is_lbl$ 

INPUT "NEU X AXIS LABEL" . X_ax i s_l bl $ 

CASE 4 

PRINT "CURRENT Y AXIS UNI TS" , Y_ax is_units$ 

INPUT "NEU Y AXIS UNITS" , Y_ax i s_uni t s$ 

CASE 5 

PRINT "CURRENT X AXIS UNI TS" . X_ax i s_uni ts$ 

INPUT "NEU X AXIS UNITS" . X_ax Is.unitsS 
CASE 6 

PRINT "CURRENT UINDOU" .Left . Right , Botton , Top 
INPUT "NEU UINDOU" .Left. Right .Botton, Top 
CASE 7 

PRINT "CURRENT X AXIS TIC SPACING" . X_t ic_sp 
INPUT "NEU X AXIS TIC SPACI NG" . X_t i c_sp 
CASE 8 

PRINT "CURRENT Y AXIS TICK SPACI NG" . Y_t ic_sp 
INPUT "NEU Y AXIS TIC SPACI NG" . Y_t ic_sp 
CASE 9 

PRINT "CURRENT TIC INTVL BTUN MAJ TICS X DIR" . X_maj _ct 
INPUT "NEU TICK INTVL BTUN MAJOR TICKS X DI R" . X_naj _ct 
CASE 10 

PRINT "CURRENT TICK INTVL BTUN MAJOR TICKS Y DIR" . Y_naj_ct 
INPUT "NEU TICK INTVL BTUN MAJOR TICKS Y DIR" , Y_naj_ct 
CASE 1 1 

PRINT "CURRENT MAJOR TICK SI ZE" . Maj_t ic 
INPUT "NEU MAJOR TICK SIZE" . Maj_t ic 
CASE 12 

PRINT "CURRENT X TICK GRID SPACING" .X_gnd_tic 
INPUT "NEU X TICK GRID SPACI NG" , X_gr id_t ic 
CASE 13 

PRINT "CURRENT Y TICK GRID SPACI NG" . Y_grid_t ic 
INPUT "NEU Y TICK GRID SPAC I NG" . Y_gr xd_t ic 
CASE 14 

PRINT "CURRENT GRID TIC INTVL BTUN MAJ TICS X DI R" . X_gr i d_na j 
INPUT "NEU GRID TIC INTVL BTUN MAJ TIC X DIR" , X_gri d_naj 
CASE 15 

PRINT "CURRENT GRID TIC INTVL BTUN MAJ TICS Y DIR" . Y_grid_naj 
INPUT "NEU GRID TIC INTVL BTUN MAJ TIC Y DIR" , Y_grid_maj 
CASE 16 

PRINT "CURRENT MAJOR TICK SIZE" ,Maj_grid_s2 
INPUT "NEU MAJOR TICK SI ZE" . Maj_gr id_s2 
CASE 17 

PRINT "CURRENT Y INCREMENT" , Y_incnt 
INPUT "NEU Y INCREMENT" .Y_incnt 
CASE 18 
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3350 

33G0 

3370 

3380 

3390 

3391 
3400 
3410 
34Z0 
3430 
3440 
3450 
34G0 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3G10 
3G20 
3G30 
3640 
3650 



PRINT "CURRENT X I NCREMENT" , X_i ncnt 
PRINT "NEU X INCREMENT" ,X_incnt 
CASE 19 

PRINT "CURRENT LINE TYPE" ,Line_nun 

DISP "NEU LINE T YPE( 1 =SOLI D , 4=DRSH , 5=BIG DASH , G=DftSH- DOT" ; 
INPUT Line_nun 
CRSE 20 

PRINT "CURRENT HEADING SUB LABEL" ,Sub_head_Ibl$ 

INPUT "NEU HEADING SUB LABEL" , Sub_head_IbI$ 

CASE 21 

PRINT "CURRENT X AXIS NUMERIC I MAGE" , X_char$ 

INPUT "NEU X AXIS NUMERIC IMAGE" . X_char$ 

CASE 22 

PRINT "CURRENT Y AXIS NUMERIC I MAGE" , Y_char$ 

INPUT "NEU Y AXIS NUMERIC I MAGE" , Y_char$ 

END SELECT 

INPUT "DO YOU UISH TO CHANGE ANOTHER PARAMETER? Y=YES",Dec5$ 

IF Dec5$="Y" THEN 
GOTO 2510 
ELSE 

GOSUB PIot_data 
END IF 
RETURN 

I 

( A A AA A A A A aA M. a. a A./ ■> A A. A. A A ■» AA - A. A.A A^ M. A AA A. A M. M. A A. ^A A ^ A A ** A -* J A. A-A. A A Ay -> A ^ A A A. A A A A A A A. A A_A_ > A, -A^ A A -A^ 4.A A. 4A aA AA -A ^ A A A A A.A 

ttnittitta tnt ti iTiriritc Pitiitrnitttiiiritatittc trttit 

xV vx yy " V ▼▼ ▼» ^ i yy tT ▼▼ ft t T ▼▼ tT " ▼ T " T ¥ ▼T 

! END THl PROGRAM 

I AA A A A A A A. A > A , j \ J j AA AA A ^ .A. A A A A A * . * A, A A. A A.A A a a a a ^ A. A A, A AA AA AA A_A A_A AA A-A AA a1 Ja AA AA AA Aa AA AA A_A A^ Aa A^ AA AA A.A Aa A A XA, AA A^ A, A A A 1 A A_a. 

I IT n if fTtrtfiTtTTT7tTTTT gtinirTtiTir n iririT irirTiixiitiiriiiTn rrTrnnnii trirn n CTTiTHTTTrTt iiiTnniiiTT ri ririiiTTXTrxf 

Tf tHT ‘TT^rr TT ^ f f ¥ • Tf yr ▼ “ T^T" f f TT II If ▼ " "T TT TT ?f Wl •* ▼? "T ^T Wt • • tt f f r'1 

! 

End_program: ! 

GRAPHICS OFF 

PRINT "END OF PROGRAM" 

END 
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